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Abstract 



We study isospin breaking instanton corrections to the operator product 
expansion of the nucleon correlation functions. After a comparison with quark 
model calculations based on the 't Hooft interaction, we examine the role of 
instantons in the corresponding QCD sum rules. Instanton contributions are 
found to be absent in the chirally even sum rule, but significant in the chirally- 
odd one. They improve the consistency of both sum rules and favor a value 
of the isovector quark condensate close to the chiral estimate. 
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I. INTRODUCTION 



Over the last years growing evidence for a significant role of QCD instantons in hadron 
structure has been collected. It originated first from models built on instanton vacuum phe- 
nomenology and recently received model independent support from cooled lattice studies 
@. Indeed, the latter show that hadron correlation functions remain almost unchanged if 
all but the instanton fields are filtered out of the equilibrated lattice configurations. 

Analytical studies of instanton contributions to the operator product expansion (OPE) 
and to QCD sum rules find a reflection of this picture in the importance of explicit instanton 
corrections in the pion j|] and nucleon Hj channels. The corrections in the nucleon channel 
show a characteristic pattern, which originates from the chirality of the quark zero mode 
states in the instanton background: they are small in the chirally-even nucleon correlator and 
in the corresponding sum rule, but significant in the chirally-odd one. Indeed, the chirally- 
odd sum rule could hardly be stabilized without instanton corrections, whereas the chirally 
even one is stable and in agreement with phenomenology even if the instanton contribution 
is neglected 0. 

An analogous pattern was found in two recent sum rule calculations of the neutron- 
proton mass difference 5Mn J^JU without instanton corrections, which also show a significant 
discrepancy between the results of the chirally even and odd sum rules. Again the former 
agrees well with phenomenology (5M^ ~ 2 MeV) , whereas the latter yields a value consistent 
with zero and thus puts the consistency of the two sum rules into question^. This analogy 
with the nucleon mass sum rules prompted us to examine instanton corrections to the isospin 
violating nucleon sum rules, which is the subject of the present paper. 



1 Attempts to reduce this discrepancy by adding a term attributed to electromagnetic corrections 
to the OPE |?J would require a substantial corresponding refinement on the phenomenological side 
of the sum rule, see ref. ||. 
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A further, closely related sum rule calculation of isospin violation in baryons without 
instanton corrections || takes a somewhat different approach. The baryon mass splittings 
are taken as input from experiment (after subtraction of the estimated electromagnetic con- 
tributions), and the two relevant isospin breaking parameters - the quark mass difference 8m 
and the difference of up- and down-quark condensates 7 - are estimated from the sum rules. 
This analysis seems to find consistency between both sum rules, at least if the difference 
5X 2 N between the neutron and proton pole strengths is fitted, and thus seems incompatible 
with the conclusions of refs. |7]J§. The fit requires, however, an unusually small value of I7I, 
about a quarter of the one estimated from chiral perturbation theory, and an uncomfortably 
large continuum contribution. We will come back to this issue below. 

The study of isospin violations in QCD nucleon sum rules can be based either on the 
nucleon correlator in an iospin-violating scalar background field M or on the difference of 
the neutron and proton correlators |7]|| . We will adopt the latter approach. In section |T| we 
calculate the leading, isospin violating instanton corrections to the nucleon correlator, and 
in section |TTT| we discuss their structure in more detail. Section |I^ contains a comparison 
with quark model calculations based on instanton-induced interactions. We point out, in 
particular, that the neglect of the vacuum sector in many of these models leads to severe 
limitations in their description of isopin violation effects. On the basis of the instanton- 
corrected nucleon correlators from section [Q| we then set up the correponding QCD sum 
rule in section and analyze it quantitatively in section VI. The final section contains a 
summary of our results and some conclusions. 

II. NUCLEON CORRELATORS 

This section describes the evaluation of small-scale instanton contributions to the nucleon 
correlators in the presence of isopin breaking. We begin with the correlation function in the 
proton channel, which is characterized by two invariant amplitudes of opposite chirality, 
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U p (q) = % J d A xe^ x (0\T Vp (x)n p (0) |0> = ^ p (q 2 ) + n 1>p (g 2 ) . (1) 

The composite operator r] p is built from QCD fields and serves as an interpolating field 
for the proton. Two such (independent) operators with minimal mass dimension (i.e. 9/2) 
can be constructed. We adopt the standard choice of Ioffe || , 

Vp( x ) = e abc [ul(x)C^u b (x) 757 M 4(^) , (2) 

(the neutron current is obtained by interchanging up and down quark fields), which allows 
for a direct comparison with the previous studies of isospin violation in nucleon sum rules 



The leading instanton contributions to the correlators can be calculated in semiclassical 
approximation, i.e. by evaluating (|1]) in the background of the instanton field and by taking 
the weighted average of the resulting expression over the quantum distribution of the instan- 
ton's collective coordinates ||. These contributions add nonperturbative corrections to the 
Wilson coefficients of the conventional OPE, with which they will be combined in section 
[V[ Isospin breaking originates in this framework from the mass difference of up and down 
quarks, 

5m = m d - m u , (3) 

and from the differences in the values of the corresponding condensates, 

_ {Q\dd-uu\0) 

7 = (dHo) ' ( ' 

The isovector quark condensate, which determines 7, is the dominant source of nonper- 
turbative isospin violation in the OPE of the correlators, since it originates from the lowest 
dimensional operators with a finite vacuum expectation value. The value of 7 has been 
estimated in a variety of approaches |]HTjJ, with results varying over almost an order of 
magnitude, — 1 x 10 -2 < 7 < —2 x 10~ 3 . The sensitivity of the baryon sum rule analysis 
to 7 can be used for an additional estimate of its value 0], which will be adapted to the 



presence of instanton corrections in section [VT]. For the quark mass difference we use the 
more accurately known standard value 5m = 3.3MeV JT0| . 



The rationale behind the semiclassical treatment of instanton contributions and the cal- 
culational strategy are analogous to those in the isosymmetric case 0, and we thus just 
sketch the essential steps here. To leading order in the product of quark masses and instan- 
ton size, instanton effects in the nucleon correlators are associated with the quark zero-modes 
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$t{*) = - ( 2+X/2 - U > (5) 

where the superscript ± corresponds to an (anti-) instanton of size p with center at Xq. The 
spin-color matrix U satisfies (a + r) U = and r = x — Xq. The zero mode contributions enter 
the calculation of the correlators through the leading term in the spectral representation of 
the quark background field propagator 

S q&V) = —^j + 0(pm q ). (6) 

The flavor dependent effective quark mass m*(p) = m q — ^7T 2 p 2 {qq) (where q stands for up 
or down quarks) in the denominator is generated by interactions with long-wavelength QCD 
vacuum fields |18|] . Quark propagation in the higher-lying continuum modes in the instanton 
background will be approximated as in |J by the free quark propagator. 

Note that both the zero and continuum mode propagators are flavor dependent. The zero 
mode part contains the effective quark mass, which depends on the current quark masses 
and on the corresponding condensates. The current quark masses enter, of course, also the 
continuum mode contributions. 

With the quark background field propagator at hand, the instanton contributions to the 
proton and neutron correlators can now be evaluated. As a first, generic result we find 
that the chirally even amplitudes U q for both proton and neutron do not receive leading 
instanton corrections. This generalizes the analogous observation in ref. || to finite current 
quark masses and condensate differences and is a consequence of using Ioffe's current. The 



isopin breaking difference of the IT^'s for neutron and proton, moreover, vanishes for all 
interpolating fields, as we will show in the next section. 

The chirally odd amplitudes ITi, on the other hand, get sizeable instanton contributions, 
and their difference for proton and neutron remains finite. For the proton correlator, to first 
order in the current quark masses and continued to euclidean space-time, we obtain 

16 t 4 n(p) 



where N c is the number of quark colors. The isoscalar part of the effective quark mass in 
the chiral limit is defined as m|$(p) = — |vr 2 p 2 (gg) with (qq)o = ({««) + (dd))/2 and the 
dimensionless ratio ( = (m u + m,j)/mj. 

The further evaluation of eq. (|?p requires an explicit expression for the instanton size 
distribution n(p) in the vacuum. Instanton liquid vacuum models [ID] and the analysis of 



cooled lattice configurations |3| have produced a consistent picture of this distribution. The 



sharply peaked, almost gaussian shape of n(p) found in ref. |T9[| can be sufficiently well 



approximated as p0 | 



n(p) = n8(p — p) , n ~ -fm 4 , p ~ -fm, (8) 

2 3 

which neglects the small half width (~ O.lfm) of the distribution. In eq. (§) we introduced 
the average instanton size p and the instanton number density n, which equals the density 
of anti-instantons. n can be approximately related |^TJ to the isoscalar quark condensate 
by the self-consistency condition n = —^m^(p)(qq)o, which quite closely reproduces the 
phenomenological value given above and allows the elimination of n in favor of the quark 
condensate. 

After performing the now trivial integration over instanton sizes, we prepare the ampli- 
tude (0) for its use in the corresponding sum rule by applying the standard Borel transform 



n(M 2 )= iim ±(Q*)r+i(-JL) n(g 2 ) (9) 



n— >oo 



' \ dQ 

(Q 2 = —q 2 ) with the squared Borel mass scale M 2 = Q 2 /n kept fixed in the limit, and obtain 



u\ n ;\M 2 



An 2 



l -^(uu)M%(M 2 p 2 ) - -Lm u p*M w I 2 (M 2 p 2 ) 
JM r lb7v z 



(10) 



in terms of the two dimensionless integrals 



h(z 2 ) = / dx , 2lA , 2 e- x2{x ~ z 2/4rl , (11) 
J Z 2/4 (x - z 2 /4) 2 

uz 2 ) = r d Xl r dx 2 x ^ z ~ 2 - x2 ^ 2 e -i(* 1+a »-^)-i (12) 



The amplitude for the neutron follows from (|10|) by interchanging up and down quark masses 
and condensates. Eq. flip]) generalizes the isospin-symmetric amplitude of ref. 0, which is 
recovered in the chiral limit. 



The two integrals (11) and (|T^) contain the instanton corrections to the Borel transformed 
Wilson coefficients of the unit operator (in I 2 ) and of (uu) (in Ii). Additional contributions 
from quark modes with momenta below the renormalization scale p of the OPE should be 
excluded from 1\ and I 2 in order to avoid double counting of the physics contained in the 
condensates. Fortunately, the instanton background field induces only one soft contribution 
up to operators of dimension 6, corresponding to a four-quark condensate in I 2 . We will 
correct for this contribution in section 0. 



III. ISOPIN VIOLATION AND INSTANTONS 

It is instructive to analyze the isospin properties and the origin of isospin breaking in the 
instanton induced amplitude (|II]) in more detail. This analysis will also lay the foundation 
for our discussion in the next section, where we clarify some crucial differences between our 
approach and quark model calculations. These differences explain, in particular, the absence 
of instanton induced contributions to 5Mn in many quark models, contrary to our findings 
in the correlator approach. 
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The gluonic sector and the quark-gluon vertex of QCD are both flavor independent. The 
structure of the interaction with the instanton background field is thus the same for up and 
down quarks and the background field propagator (|6]) is diagonal in isopin space. Its only 
flavor dependence enters through the current quark masses and condensates in m*. This has 
characteristic consequences for the instanton-induced interactions between quarks, which 
generate the correlator amplitude (|10D and can be extracted from the calculation in section 







(0| T q A , a ,a{xi) qB,p,b{yi) qcriA^) QDfrdfa) |0) 



n(p) 



(5 AB 5 CD - 5 AD 5 CB ) J d P {m ^ p) (2^ P r 

x J d 4 a;oC(r 1 )C(n 1 )C(r 2 )C(n 2 ) £ « (P L /RPab)af)(PL/RPcdU » 



SU(3) C 



(13) 



L/R 



In this expression, capital latin, greek and small latin indices refer to isospin, Dirac spin 
and color, respectively. The brackets << ... >> indicate the average of the instanton's 
color orientation over the Haar measure of SU(3) C . The chiral projection operators are 



P 



L/R 



[1 ±75)72, the distances from the instanton center are denoted 



X j 



u i — yiUi ~~ x o) 2 i the nonlocality of the vertex is contained in the functions 

( ) ^ r 2 _|_ p2^j 

and its color structure is given by the spin-color tensor 



Xq) 2 and 



(14) 
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(15) 



ab 



which makes the embedding into an SU(2) subgroup of SU(3) C explicit and contains the 
characteristic spin-color coupling. Lorentz and color covariance become manifest only after 
averaging over the color group. 

The nonlocal four-quark vertex flT3"D originates from the quark zero modes in the instan- 
ton field and was first derived by 't Hooft [0 for m* = m q . It has, separately for both quark 



chiralities, a well-known determinantal flavor structure and is thus SU(2)l X SU(2)r and in 



particular isospin symmetric?]. Even if the flavor dependent quark masses and condensates 
enter the vertex explicitly, isospin violation can thus not originate solely from the instanton 
induced interaction. In the OPE, however, this vertex generates nonperturbative contribu- 
tions to the Wilson coefficients, which can either become themselves isopin violating due to 
the finite quark mass difference or multiply isopin violating operators. Examples for both of 
these cases were found in section |I[ 

At short distances, multi-instanton contributions to the nucleon correlators are sup- 
pressed since x/R <C 1. Additional corrections, which originate from only one valence 
quark propagating in a zero mode state and generate contributions to the quark self-energy, 
are subleading in m*p and not considered in this paper. 

The structure of the instanton-induced four-point function ( O ) explains some qualitative 
features of the correlator flTUp. The definite chirality of the quark legs, inherited from the 
zero mode states, lets this vertex act only between quarks which are coupled to spin 
in the interpolating fields, and the Dirac structure of the instanton contribution to the 
nucleon correlator is thus determined by the remaining "valence" quark line. Only the 
contribution to the Wilson coefficient of the unit operator, in which this line is in a non- 
zero-mode state (approximated by the free quark propagator), can thus generate a chirally 
even amplitude. According to our discussion above, however, such contributions are isospin 
conserving and this explains why the difference of the chirally-even amplitudes of the neutron 
and proton correlators is not corrected by instantons for any choice of the interpolating field. 
The instanton contributions to the chirally even amplitudes of neutron and proton vanish 
individually only for the Ioffe current, as already pointed out. 

Another characteristic feature of the correlator ([H]) follows directly from the flavor struc- 
ture of the 't Hooft vertex. Since only one quark pair of each flavor can take part in the 



2 In addition it breaks the axial ?7a(1) symmetry, which is a reflection of the Adler-Bell-Jackiw 
anomaly and a celebrated instanton effect 17]. 
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zero mode induced interaction ([13]), the valence quark in the non-zero-mode state in the 
proton (neutron) correlator must be an up (down) quark. This explains why we find only 
contributions proportional to the up quark mass and condensate in eq. (|10|). 

To summarize, instanton induced interactions contribute to isopin breaking in the oper- 
ator product expansion of the nucleon correlators (in the framework of our approximations) 
in two distinct ways: they correct the Wilson coefficient of the unit operator, which becomes 
isopin dependent due to the difference of the current quark masses, and they contribute to 
the coefficients of the isospin violating operators uu and dd. Both of these corrections affect 
only the chirally odd amplitude of the correlators. 



IV. COMPARISON WITH QUARK MODELS 

Instanton induced interactions have been included in several quark model calculations of 
mass splittings in baryon iso-multiplets. It is useful to compare the results of such calcula- 
tions to those of our approach. To be specific, we will base this comparison on studies in 
the MIT bag model p3"| , p3| , which deals with relativistic, light quarks and is in this respect 
similar to the correlator approach. Most of our conclusions, though, will apply to a wider 
range of quark models. 



Bag model calculations of hadron mass shifts due to instantons |25[] are based on a 
localized version of the 't Hooft interaction, cast into the form of an effective lagrangian. 
Indeed, the pointlike limit of the vertex ( |13"D (in Minkowski space) is reproduced by the 
lagrangian [TF 



4 n { 

3 {m* u p){m* d p) f^ R { 



3 

F 32 



3 

[u R X a u L )(d R X a d L ) - -{u R X a o^u L ){d R X a a^ u d L ) 



(16) 



which is obtained from flT3| ) by amputating the external quark propagators, neglecting the 
nonlocality due to the finite instanton size, performing the average over the color orientation 
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of the instanton]], specifying the instanton density in the form @ and continuing back to 
Minkowski space-time. 

The shrinking of the instanton vertex to its pointlike limit will probably not cause sig- 
nificant errors in bag model results for low-lying hadrons. This is because quarks in the bag 
can separate up to the diameter 2R ~ 2fm, so that their inverse momenta in the ground 
state are considerably larger than the average instanton size which characterizes the extent 
of the vertex. Note, however, that the nucleon correlator in QCD sum rule calculations is 
probed at an order of magnitude smaller distances, where the details of the short distance 
dynamics and thus the nonlocality of the vertex become important. 

Besides the structure of the instanton induced quark interaction, bag calculations share 
some other common features with the nucleon correlators ([]]), notably in the construction of 
the nucleon states. The spin, color and flavor structure of the bag model (i.e. SU(6)) proton 
state, 

|P> T) = 77f| e ^ c [( M a|4 T - uJ^H^ct] |0> (17) 



(Arrows indicate the value of the total spin projection j% of the quarks (in the bag ground 
state) and of the proton.) and the corresponding one for the neutron (which is obtained from 
— \p) by interchanging up and down quarks) is essentially identical to that of the interpolators 
(0). This is of course just a consequence of the fact that both are constructed to carry nucleon 
quantum numbers, which ensures that they have identical properties under Lorentz, color, 
isospin and the standard discrete transformations. 

Despite these similarities, bag model calculations do not find any instanton contribu- 



3 As long as this lagrangian is evaluated only in color singlet states, one could of course skip the 
color averaging and use the neither Lorentz nor SU(3) C invariant version instead, with identical 
results. 
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tion to the proton-neutron mass difference^ p3| , |24"|| . The neutron and proton mass shifts 
induced by (|T6|) are evaluated in first order perturbation theory between the SU(6) states. 
A straightforward calculation gives 

( P \ - [ d*xC inst {x) \P) = & 2 P 3 ) 2 ~ ( v4( ^\ (^) 2 
J 3 {m* u p){m* d p) V 8tt J 

x J d 3 x {5 {a\ - alf u \a\ - a 2 ) (d) + 4 (a? + a 2 )W(a 2 + a 2 2 )^ + 20 a^a^a^af) , (18) 

(the integrations extend over the bag volume) where the ground state quark wavefunction 
in the bag is written as 




wi th a[ q) = y^Sj'o (¥) and a { 2 q) = Ji (^). (-R is the bag radius, N q a nor- 

malization constant, \n are the Pauli spinors and u q , K q are the energy and (dimensionless) 
momentum quantum numbers of the quark ground state.) The result for the neutron matrix 
element can be immediately infered from (ji~8"D by exchanging up and down quark operators 
in the proton states. Since £i ns t is symmetric under u <-» d, one can as well exchange u 
and d everywhere in the expression for the proton matrix element. Eq. fll8|), however, is 
manifestly invariant under this exchange and therefore the difference of the matrix elements 
indeed vanishes: 

SMpfcg = (p| J d 3 x £ mst \p) - (n\ J d 3 x £ inst \n) = 0. (20) 

This is in contrast to the result for the Borel transformed nucleon correlators, 

n ^f (M 2 ) - n^'(M 2 ) ^ 0, (21) 

which can be translated via a QCD sum rule into a finite instanton contribution to the mass 
difference (see below). 



4 as long as one-zero-mode corrections (see above) are neglected 
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Of course, the bag model contains ad hoc assumptions on quark confinement, breaks 
chiral symmetry explicitly and differs also in other aspects from the model-independent 
correlator approach. One would thus not expect the results to agree quantitatively. It is at 
first surprising, however, to find a qualitative difference, namely the exact absence of any 
instanton mediated mass shift 5M^f ag . In view of the similarity of the interactions fll3|) and 
(|16|) one is lead to search for the origin of this difference in the description of the nucleon 
states. And indeed, here the two approaches differ crucially. 

A first difference is that the quarks in the bag matrix elements are restricted to their 
ground states with total spin j = 1/2 (i.e. to / = or 1), whereas the interpolators can create 
quarks in all orbital angular momentum states. In fact, the created pointlike wavepacket has 
overlap with the whole tower of excited states carrying nucleon quantum numbers, including 
states in the many-particle continuum. Experience from QCD sum rules || and from lattice 
data [^,0 shows, however, that already at rather small distances a main contribution 
to the nucleon correlator comes from the nucleon ground state. This tendency is further 
enhanced by the Borel transform, which exponentially suppresses contributions from higher 
lying states. One does therefore not expect these states to contribute significantly to a finite 
value of the difference ( pTj) in the fiducial Borel mass domain (see below), let alone to be its 
only cause. 

Indeed, the crucial difference between the bag and correlator results can rather be traced 
to the description of the nucleon ground states themselves, and in particular to their flavor 
content. While the SU(6) states ( |17|) (as well as the interpolating fields (§)) have good 
isospin, this is not the case for the states i] Pjn \0) which are created by the interpolators and 
studied in the correlator ([!]). 

Virtual "sea" quarks and other perturbative and nonperturbative vacuum fluctuations 
give these states a much richer (and more realistic) flavor structure. They inherit, in partic- 
ular, isospin breaking components from the vacuum fields. The short distance part of this 
nontrivial flavor content is captured both in the OPE and in the instanton corrections and 
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causes ( pT]) to be finite. At very short distances it originates from the quark mass differences, 
and at larger distances it enters predominantly through the difference between up and down 
quark condensates. The neglect of both of these ingredients in the isospin structure of the 
bag wavefunctions leads, on the other hand, to the symmetry of the proton matrix element 
([Of) under exchange of the two quark flavors in the states ( |17D and thus to the absence of 
an instanton-induced mass difference, eq. (|20|). 



In ref. [^4| an attempt was made to include long-wavelength vacuum fields into the bag in- 
terior. Long and short distance physics inside the bag, however, cannot be reliably separated. 
In particular, such a scale separation (which is indispensable to control the interactions with 
the vacuum fields) cannot be based on a short distance expansion, which would in fact badly 
diverge at distances of the order of the bag radius R ~ ^qcd- Higher order interactions 
with the background fields (leading to contributions from higher dimensional condensates) 
are thus not suppressed, and it is hard to see how their neglect can be justified and how 
double counting of quark physics can be avoided. In addition, both the contributions of 
the interactions with the long-wavelength background fields and with the instanton to the 
matrix elements are calculated independently to leading order and then added. Combined 
effects of instantons and the other vacuum fields, as described by the OPE, are therefore 
still lacking^ and eq. (|20|) remains to hold. 

The large distance scales over which quarks can interact are a generic bag model problem, 
since clearly not all nonperturbative physics can be absorbed into the boundary conditions. 
For QCD sum rule calculations, on the other hand, a reliable description of the correlation 
functions up to distances of about 0.2 fm is usually sufficient. At these rather small distances 
the long- wavelength physics can still be controlled in a model independent way by only a 
few generic and physical parameters, the low-dimensional condensates. 

Instanton physics supplies, in fact, yet another example for the problems with treating 



-"except for the factor (m^m^) 1 in the instanton-induced interaction (16) 
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quark interactions in the large bag interior. The neglect of multi-instanton effects (be- 
yond the mean field level) is a basic assumption underlying the interactions (|T3|) and (fl6|). 
This approximation can be justified in the correlator for distances x -C R, but hardly at 
the scales set by the bag diameter. Instanton liquid simulations [28j indeed confirm that 
multi-instanton effects become important at distances of the order of the average separation 
between instantons, x > R ~ 1 fm. 

To conclude, quark and bag model calculations which restrict the evaluation of instanton- 
induced baryon mass shifts to calculating the expectation value of the effective 't Hooft 
lagrangian between SU(6) eigenstates miss important sources of isospin asymmetry, in par- 
ticular from isospin violating vacuum fields. This physics, which significantly affects the 
estimates of mass splittings in baryon iso-multiplets, is however captured in the instanton- 
corrected OPE of the nucleon correlators. 



V. ISOPIN VIOLATING NUCLEON SUM RULES 

In this section we combine the instanton contributions to the nucleon correlation functions 
with the conventional operator product expansion [[/J and set up the QCD sum rule for the 
difference of the neutron and proton correlators. 

Since the average instanton size p is smaller than the inverse renormalization point 
~ 0.4 fm of the OPE, the major part of the instanton corrections will contribute to 
the Wilson coefficients. These corrections can be directly added to the standard OPE, since 
they originate from nonperturbative physics which was not previously accounted for. 

The integral I2 in eq. (|T0|), however, contains besides the instanton contribution to the 
Wilson coefficient of the unit operator also a soft part, as pointed out in section |TJ. It 
originates from the region in loop momentum space where the hard external momentum Q 2 
is carried exclusively by the quark line which is not participating in the zero-mode induced 
interaction. No such contribution is contained in Ji, since here the third quark line interacts 
with the quark condensate and thus does not carry momentum. 
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The soft part of I 2 represents an instanton contribution to the four-quark condensate. 
Fortunately this is the only operator up to dimension 6 which receives such a correction. 



Indeed, a general theorem |29| severely limits the number of condensates in the OPE of 
hadronic correlators which can be induced by self-dual background fields. 

In order to prevent double counting of long-wavelength physics, the four-quark condensate 
terms in the OPE and the instanton contributions have to be adapted before combining them 
with the other OPE terms of III. As in ref. || we will neglect the comparatively small OPE 
contribution to (uudd) and keep instead the contribution induced by the 't Hooft vertex in 
the limit of vanishing external momenta. A more accurate procedure, namely to subtract 
explicitly the part of the instanton contribution which originates from momenta below the 
renormalization scale, will be described elsewhere Four-quark condensates of the type 



(uuuu) and (dddd), on the other hand, do not receive single-instanton contributions and 
remain unchanged. 

After implementing the above modification, we can combine the instanton part (0) with 
the OPE of the chirally odd sum rule of ref. 0. Taking the difference of the neutron and 
proton sum rules and transfering the continuum contributions to the left-hand (i.e. OPE) 
side, we obtain 
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M' 

■2\- X j±<S.\l x - X 2 N 5M N - 5X 2 N M N 



M%/M* 



1 



■(qq)oSie-^ M2 5 Sl , (22) 



4tt 2 

where Mn = (M p + M n )/2 and Xn = (A p + A n )/2 denote the isoscalar nucleon mass and 
coupling to the interpolating field. (We neglect the small gluon condensate contribution.) 
The isospin violating differences of the overlap and threshold parameters are 5\ 2 N = — X 2 
5si = s Vn — Sip, and the factor L -8 / 9 , with L = ln(M 2 / Aq CB ) / ln(fi 2 / Aq CB ) and A QCD = 
150 MeV, accounts for the anomalous dimensions of the composite operators and sets their 
renormalization point to /i = 500 MeV. 

The contributions from the continuum, starting at the effective threshold Si, are as 

16 



usual combined with the leading OPE term and described by the functions E\ = 1 — 
e~ Sl/M2 (jfe + l) and E 2 = 1 - e - s ^ M ' 2 + + l) @. Their definitions are identi- 
cal to those in the individual sum rules. Note that additional terms in the difference sum 
rules, proportional to 5si and 5s q (see below), originate from the continuum terms of the 
individual sum rules. They do not correspond to the cut structure (i.e. to the leading OPE 
behaviour) of the difference sum rule, however, and are thus not needed to match the large-s 
behavior of the OPE. 

With the standard definitions a = —4iT 2 (qq) , and X 2 N = 327r 4 A^, our sum rule eq. ( |22|) 
now assumes its final form 



M 8 6mE 2 L~ 8 / g + M 6 7 a£i + pmM 2 a 2 - M^al^z 2 ) - ^5mp 4 M 12 1 2 (z 2 )L- 8/g 



= \ 2 N M 2 N 5M N - \^f5M N + ^M N j M 2 + as^e^-^^ds, . (23) 

The corresponding sum rule from the ^ structure |7|-||] is unaffected by leading instanton 
corrections, 



e M%/M* 



-aM 4 SmE L- 4 / 9 - -M 2 ia 2 L^ + ^5mM 2 aL~ 8 / 9 + ^a 2 ^ 2 ' 27 

3 6 3 

= \ 2 N M N 6M N - ^M 2 + \(s* + ^ M 2 e-^ M2 ^ M2 L- A ' 9 5s q . (24) 

The two parameters b = (glG^G 1 ^) = 0.5 GeV 4 and m 2 , = (g s qa ■ Gq)/(qq) = 0.8 GeV 2 are 
fixed at their standard values and E = 1 — t~ Sq l M2 . 

Above we have written both sum rules, eqs. (|23D and (0), i n their most general form, 
which allows for independent values of the isosymmetric continuum thresholds Si and s q . 
Below we will, however, follow the standard practice and set si = s q = sq. 

At this point it might be useful to recall the main assumptions and approximations which 
went into the OPE of these sum rules and into the parametrization of their phenomenological 
sides. (For more details see [0,^]-) The short-distance expansion is carried out up to 
operators of dimension six and the perturbative part of the Wilson coefficients is calculated 
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to leading order in the strong coupling a s . Less systematic uncertainties arise from the 
not precisely known values of the condensates and from the standard factorization of the 
four-quark condensates. 

As is common practice in sum rule calculations, the Wilson coefficients are calculated 
without explicit infrared cutoff since at scales up to about 0.5 — 1 GeV non-perturbative 
contributions strongly dominate over the perturbative ones. The explicit removal of the 
latter becomes therefore practically unnecessary f3lfl . For the same reason, the condensates 
(and the quark masses) depend in the above range rather weakly on the renormalization 
scale. Even without explicitly specifying the infrared regularization scheme of the Wilson 
coefficients scale-dependent quantities are understood to be taken at a fi in the above range, 
and we use specifically fi = 0.5 GeV. 

On the phenomenological side the main assumption is that of local duality between the 
hadron and quark-gluon descriptions of the continuum. It has been found to work well in 



many sum rule studies and also in recent lattice simulations of point-to-point correlators [26 



In our context it is put to a harder test since we consider differences of two spectral functions. 
Here even more than in the single nucleon sum rules the exponential Borel suppression of 
continuum contributions is important in order to increase the sensitivity of the sum rules to 
the ground state contributions. 

In the numerical evaluation of the sum rules the upper limit of the Borel interval is 
determined such that contributions from the continuum do not exceed a given percentage 
of the full OPE contribution. Otherwise the sum rules would be relatively less sensitive to 
the pole contribution of interest and a good fit quality would become a trivial consequence 
of continuum domination (instead of being a consistency criterion), since the continuum is 
modeled after the leading OPE behavior. Moderate continuum contributions are therefore 
a necessary condition for reliable sum rules, and in the next section we will check how these 
contributions are affected by the instanton terms. 
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VI. QUANTITATIVE SUM RULE ANALYSIS 



The quantitative analysis of isospin violation in the nucleon sum rules aims at determin- 
ing the isospin breaking parameters on the phenomenological side from the best fit to the 
"theoretical" left hand side. Taking all the other parameters from the standard, isosym- 
metric nucleon sum rules or from experiment, it would still require a four-parameter fit to 
determine SM^, SX%, Ss q and 8s\ independently. Limitations in the parametrization of the 
spectral densities and approximations on the theoretical side would, however, make such a 
fit unstable both with and without instantons. 

In order to reduce the number of fit parameters, one is thus led to either fix the only 
phenomenologically known one, SM^, at its experimental value M or to make assumptions 
relating at least two of the remaining isospin breaking parameters. The authors of ref. [[7], 
for example, assume 5s q = 5s± in their analysis or, alternatively, neglect differences in the 
effective continuum of proton and neutron channels entirely, i.e. Ss q = 5si = 0. Since such 
assumptions lack theoretical foundation, the associated errors can not be reliably estimated 
or controlled. We thus prefer to follow the approach of ref. @, taking 5M™ n ~ elm = 2.05 ± 
0.30 MeV as input from phenomenology. This value is derived from the experimental mass 
difference dM^f 1 = 1.29 MeV |HJ by subtracting the electromagnetic contribution 5M^ m = 
-0.76 ±0.30 MeV ||. 

For the isoscalar nucleon mass and quark condensate we use the standard values = 
940 MeV and (qq)o = —(225 MeV) 3 . The residuum of the isosymmetric nucleon pole, = 
1.8 GeV 6 , and the isospin average of the continuum threshold, so = 2.2 GeV 2 , are obtained 
from the instanton-corrected nucleon mass sum rules M in the same Borel window as the 
one used below. 

The isospin-breaking parameters SX 2 ^, Ssi and 5s q are then calculated by minimizing the 



difference between the left- and right-hand sides of the sum rules (|23|) and fl24j) under the 
logarithmic measure 5 of ref. || in the Borel-mass region 0.8 GeV 2 < M 2 < 1.4 GeV 2 . 
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We performed this minimization for various values of 7 in the range —1 x 10~ 2 < 7 < 
—2 x 10~ 3 discussed in section We find a better agreement between both sum rules towards 
larger (and more conventional) values of I7I in this interval, whereas the analogous study of 
||, which neglects instanton contributions, prefers an unusually small value, I7I = 2 x 10~ 3 . 
For 7 = — 1 x 10~ 2 , in particular, our best fit between LHS and RHS results in the parameter 
values 

6Xjf = -2.1 x 1(T 4 GeV 6 , 8 Sl = -1.7 x 1(T 2 GeV 2 , 5s q = 1.03 x 10~ 3 GeV 2 . (25) 

The Figures 1-3 show different aspects of this fit. In order to compare the fits of the 
optimized U q and Hi sum rules, we transfer all but the first two terms on the RHS of the 
rig sum rule (^) to the left and we rewrite the sum rule (^) analogously, so that the 



~ (5A 2 

same two terms, \ 2 n Mn5Mn ^"M 2 , remain on its RHS. Figure 1 compares this RHS 

(continuous line) with the modified LHS of the n q (dotted line) and Hi (dashed line) sum 
rules. 

In Figure 2 we plot the resulting neutron-proton mass difference SM^ as a function of 
M 2 , obtained by solving both optimized sum rules for 8M^(M 2 ). These curves show an 
extended stability plateau, which confirms the satisfactory agreement between the two sum 
rules. Indeed, this Borel mass independence of observables is the only intrinsic consistency 
criterion for the sum rules. 

In order to compare the relative size and behavior of the OPE and instanton contributions 
to the Hi sum rule (|23f) (recall that the U q sum rule does not receive instanton corrections), 
we display both of them separately, as well as their sum and the fit to the optimized RHS of 



eq. (|23[) , in Figure 3. The instanton terms reach almost the magnitude of the perturbative 
and power terms and play clearly an important role in determining the sum rule results. The 
usual practice to neglect these contributions seems thus unjustified. 

It is also instructive to compare our results of Figures 1 and 2 with the analogous curves, 
but calculated without instanton corrections. Recall that in this case the (uudd) part of 
the four-quark condensate has to be restored in eq. (p3|), which changes the factor 4/3 in 
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its Wilson coefficient to —2/3. As already noted a smaller absolute value of 7 is favored in 
this case, and the curves in Figures 4 and 5 were obtained by optimizing the sum rules with 
7 = —2 x 10~ 3 . Up to small corrections from the neglected 8-dimensional condensates, they 
correspond to the onesf] analyzed in ||. 

From Figure 4 it is also clear that rather different values of the isospin breaking parame- 
ters (5X% = 1.4 x l(T 2 GeV 6 , Ss q = 7.0 x l(T 3 GeV 2 , 5s 1 = 1.2 x 10" 2 GeV 2 ) are required to fit 
phenomenological and theoretical sides as long as instanton contributions are neglected. The 
difference between the pole strength of neutron and proton, in particular, becomes about 
two orders of magnitude larger and changes sign. 

More importantly, however, the small modulus of 7 preferred by this fit has an unwelcome 
consequence. Closer inspection of the sum rules reveals that decreasing values of I7I lead 
to increasing contributions from the continuum relative to the power corrections in the 
optimized sum rules. Indeed, the parameter values used above correspond to a continuum 
contribution of 90% in the chirally-odd sum rule (and about 37% in the chirally-even sum 
rule). This continuum domination casts serious doubts on the reliability of the chirally-odd 
sum rule, even if fit quality and stability seem satisfactory (cf. Fig. 5). In both instanton- 
corrected sum rules, on the other hand, the continuum contributions remain moderate (about 
20%). 



It is interesting to note that another recent sum rule analysis of 7 [34], which is based on 
the mass splitting in the D and D* isospin dublets, also finds a small value I7I ~ 2.5 x 10~ 3 , 
which is close to the result of Ref. ||. Since it is derived from an independent sum rule the 
discrepancy with the result of chiral perturbation theory might have a different origin in the 
D meson channel. This issue and the role of instanton corrections in this channel deserve 
further investigation, which will the be subject of a forthcoming publication | |35| |. 



6 The sum rule of ref. H contains an error in the coefficient of the four-quark condensate which 
we have corrected. 
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VII. SUMMARY AND CONCLUSIONS 



In this paper we study the role of instantons in the dynamics of isospin violation, as it 
manifests itself in the short distance expansion of the nucleon correlation functions. Isospin 
breaking effects lead to differences between the neutron and proton correlators, which can 
be translated via dispersion relations into isopin violating vacuum and nucleon parameters. 

The isospin-breaking instanton corrections to the nucleon correlators show several char- 
acteristic qualitative features. As a consequence of using Ioffe's interpolating field, instanton 
contributions are absent in the chirally-even amplitudes. Moreover, the difference of these 
amplitudes for neutron and proton is not affected by instanton corrections for any choice of 
the interpolating field. 

The chirally odd amplitude, on the other hand, receives instanton contributions of almost 
the magnitude of the standard OPE terms, as in the isosymmetric case. They correct the 
Wilson coefficients of the unit operator and of the quark condensates. The difference between 
the neutron and proton amplitudes is, in fact, mainly generated by the quark condensate 
terms, i.e. by isopin violating quark modes in the vacuum. 

This confirms the general expectation that isospin breaking effects in hadrons are phys- 
ically subtle not only because they are small, but in particular because they depend sensi- 
tively on non-valence-quark physics. This is a challenging and little tested regime for hadron 
models, which often neglect vacuum effects alltogether and thus miss important sources of 
isospin asymmetry. Bag (and other quark) model calculations which evaluate instanton- 
induced quark interactions between SU(6) states with good isospin fail, for example, to 
find instanton contributions to the neutron-proton mass difference. The instanton-corrected 
OPE, on the other hand, contains vacuum physics at short distances and thus provides a 
more reliable and model-independent basis for the study of isospin breaking effects. 

The link between the correlators and nucleon properties is established by dispersion 
relations and takes the form of two QCD sum rules for the difference of the neutron and 
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proton amplitudes. In adopting an approach for their quantitative analysis one has to decide 
between several alternatives. Taking the RHS to be the difference of the conventional pole- 
continuum ansatze for the neutron and proton, it contains four isospin-breaking parameters 
which cannot be determined independently from a stable fit, even if instanton corrections are 
taken into account. In this situation one can either assume relations between these parame- 
ters or one can fix the only phenomenologically known one, the nucleon mass difference, at 
its experimental value. We adopt the latter approach since it does not introduce additional 
assumptions with uncontrolled theoretical errors. 

The resulting sum rules, including the instanton corrections, are stable and receive only 
moderate (~ 20%) continuum contributions. This is a clear improvement over the analogous 
analysis without the instanton terms, where the continuum dominates. At the same time, the 
instanton contributions reduce the difference between the nucleon pole strengths and enhance 
the corresponding shift in the effective continuum thresholds. Moreover, and perhaps most 
importantly, the optimization of the sum rules with direct instanton effects favors larger and 
more standard values for the modulus of the isovector quark condensate, |7| ~ 1CT 2 , which 
are close to those found in the chiral analysis. 

We also tested an alternative approach towards the sum rule analysis. In this case 
the differences between proton and neutron continuum thresholds in both sum rules were 
assumed to be equal and the neutron-proton mass difference 5M N was determined from the 
fit. Inclusion of the instanton part allows a consistent fit of both sum rules, which seems 
otherwise impossible. The value of 5Mn is then, however, overestimated by about 80 %. 
This puts the initial assumption of an equal deviation of neutron and proton thresholds from 
the isoscalar position into question. It also supports our preference for the analysis method 
discussed above, which does not require ad-hoc assumptions to relate fit parameters. 

We conclude that instanton corrections play a significant role in the analysis of isospin 
breaking in nucleon sum rules. They, for example, strongly affect the results for the difference 
between the nucleon pole strengths and for the shifts in the effective continuum thresholds. 



23 



In addition, the instanton corrections enhance the internal consistency of the sum rules and 
predict a larger and more standard value for the modulus of the quark condensate difference, 
|7| ~ (0.8 - 1) x 1CT 2 . 
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FIG. 1. Best fit of the RHS (continuous line) of the sum rules to the the LHS of the Tl q (dotted 
line) and 111 (dashed line) sum rules. 



FIG. 2. The neutron-proton mass difference as a function of the Borel mass from the optimized 
H q (continuous line) and ITi (dashed line) sum rules. 



FIG. 3. Instanton and OPE contributions to the LHS of the IIi sum rule. Their sum (dashed 
line) is fitted to the RHS (continuous line). 



FIG. 4. Same as Fig.l for the sum rules optimized without instanton contributions and with 
7 = -2 x 1(T 3 . 



FIG. 5. Same as Fig.2 for the sum rules optimized without instanton contributions and with 
7 = -2 x 10~ 3 . 



28 



o 

bo 



o 

CO 



£\3 

< 

ro 



CO 



cn 



ro 
b 



CO 



CO 

b 



CO 



b 



cn 



cn 
b 




o 

bo 



o 

CO 



o 
b 



£\3 

< 

ro 



CO 



M n -M p (MeV) 

ro 

01 b en b 



ro 
cn 



GO 

b 




